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Several  large  research  programs  focused  on  the  Adriatic  Sea  in  winter  2002- 2003,  making  it  an  cs- 
,  citing  place  for  sediment  dynamics  modelers  (figure  I  i.  Investigations  of  atmospheric  forcing  and 
oceanic  response  ( including  wave  generation  and  propagation,  water-mass  formation,  stratifica¬ 
tion,  and  circulation),  suspended  material,  bottom  boundary  layer  dynamics,  bottom  sediment,  and 
small-scale  stratigraphy  were  performed  by  European  and  North  American  researchers  participating 
in  several  projects.  1  he  goal  ol  l.uroS  FRAI AIORM  researchers  is  to  improve  our  ability  to  under 
stand  and  simulate  the  physical  processes  that  deliver  sediment  to  the  marine  environment  an, I  gen¬ 
erate  stratigraphic  signatures.  Scientists  involved  in  the  Po  and  Apcnnine  Sediment  Transport  and 
Accumulation  !  PAS  I A !  experiment  benefited  from  other  major  research  programs  including  ACE 
(Adriatic  Circulation  Experiment),  IXM.CE  VITA  •:  Dynamics  of  Localized  Currents  and  Eddy  Vari¬ 
ability  in  the  Adriatic),  LACE  ( the  Croatian  Last  Adriatic  Circulation  Experiment  project !.  WISE 
(West  Istria  Experiment  i, and  ADRICOSM  (Italian  nowcasting  and  forecasting)  studies. 
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CHALLENGING  DEPOSITIONAL 
PATTERNS 

A  motivation  for  studying  the  Adriatic 
Sea  is  the  shape  of  Holocene  deposits 
and  the  separation  between  the  sources 
of  sediment  supply  and  regions  of  long¬ 
term  sediment  deposition  (see  Figure 
1  in  Nittrouer  et  al.,  this  issue).  The  Po 
River  is  the  largest  single  source  of  sedi¬ 
ment  to  the  Adriatic  and  provides  about 
one  quarter  ( 12  Mt/yr)  of  the  ~  47  Mt/yr 
sediment  delivered  to  the  Adriatic  Sea. 
Several  short,  steep  rivers  draining  the 
rapidly  uplifting  Apennine  mountains 
have  the  potential  for  delivering  high 
sediment  loads  at  concentrations  that 


make  the  river  waters  more  dense  than 
seawater  (hyperpycnal  flows).  The  Apen¬ 
nine  rivers  contribute  about  32  Mt/yr, 
but  since  World  War  II  these  rivers  have 
been  highly  regulated,  and  now  are  less 
likely  to  deliver  sediment  in  concentra¬ 
tions  high  enough  to  trigger  hyperpy¬ 
cnal  flow.  Geophysical  surveys  indicate 
that  about  180  km3  of  sediment  have 
accumulated  along  the  western  margin 
of  the  Adriatic  in  elongate  depositional 
features  termed  clinoforms,  with  the 
locus  of  deposition  in  40  to  50  m  water 
depths  (Cattaneo  et  al.,  2003;  Cattaneo 
et  al.,  this  issue).  Local  accumulations 
occur  near  sediment  sources  at  the  Po 
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delta  and  along  the  Apennine  margin, 
where  they  thicken  southward,  presum¬ 
ably  reflecting  the  cumulative  contri¬ 
bution  of  multiple  rivers.  However,  the 
thickest  accumulations  are  located  north 
of  the  Gargano  promontory,  nearly  400 
km  from  the  Po  and  100  km  from  other 
significant  sources  of  sediment.  How  can 
numerical  models  of  waves,  circulation, 
bottom  boundary  layer  dynamics,  and 
sediment  transport  be  used  to  under¬ 
stand  these  patterns? 

OLD  AND  NEW  TRANSPORT 
MECHANISMS 

Traditional  views  of  sediment  transport 
on  continental  shelves  have  focused  on 
wave-induced  resuspension  and  current- 
induced  suspended-sediment  transport. 
In  many  of  the  well-studied  systems,  the 
combination  results  in  dispersal  systems 
that  combine  advection  and  diffusion  in 
varying  amounts  (Swift  et  al.,  1972).  In 
most  cases,  however,  the  locus  of  depo¬ 
sition  is  close  to  the  source  of  sediment 
supply,  and  accumulation  rates  decrease 
with  distance  in  a  pattern  that  reflects 
the  long-term  correlation  among  sedi¬ 
ment  concentrations  and  current  ve¬ 
locities.  Observations  made  during  the 
STRATAFORM  study  of  the  Eel  River 
margin  offshore  California  revealed  a 
paradox  that  challenged  these  traditional 
views.  There,  measurements  during  the 
fifth  largest  flood  in  93  years  indicated 
that  rapid  deposition  of  flood  sediments 
occurred  beneath  the  river  plume  in 
shallow  water,  but  long-term  accumula¬ 
tion  was  centered  in  deeper  water  (Som- 
merfield  and  Nittrouer,  1999;  Wheatcroft 
and  Borgeld,  2000).  Wave-induced  remo¬ 
bilization  of  the  muddy  flood  deposits 
and  subsequent  downslope  density  flows 
were  key  mechanisms  of  cross-shelf  sedi- 
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Figure  2.  (a)  Mud  initially  deposited  in  shallow  water  (20-40  m)  beneath  the  Eel  River  plume  is  remobilized  by  waves 
and  moved  downslope  as  a  gravity  current,  producing  a  final  deposit  in  deeper  water  (50-60  m),  accord  ing  to  these 
model  simulations.  The  first  field  observations  of  this  process  were  made  during  the  STRATAFORM  project,  and  fluid 
muds  were  again  observed  on  the  flanks  of  the  Po  River  delta  during  EuroSTRATAFORM.  The  model  (b)  indicates  that 
these  processes  occur  in  shallower  depths  in  the  Adriatic  Sea,  where  waves  are  smaller.  Note  the  deposits  are  plotted 
with  ~10x  vertical  exaggeration.  Nittrouer  et  al.  (this  issue)  show  maps  of  the  final  location  and  thickness  of  fluid  mud 
deposits  after  the  Po  flood  of  2000  predicted  from  a  two-dimensional  model. 


ment  transport  (Traykovski  et  al.,  2000) 
(Figure  2).  However,  the  Eel  margin  has 
large  waves  and  high  rates  of  sediment 
supply.  Another  important  objective  of 
the  EuroSTRATAFORM  program  was  to 
determine  whether  wave-induced  den¬ 
sity  flows  occur  in  more  moderate  envi¬ 
ronments  such  as  the  Adriatic  Sea. 

SCALING  FROM  EVENTS 
IN  A  SINGLE  SEASON  TO 
GEOLOGICAL  TIME  SCALES 

One  of  the  big  challenges  in  sediment 
modeling  is  scaling.  We  do  not  have  the 
computer  power  to  resolve  processes  at 
the  finest  scales  and  also  cover  the  nec¬ 
essary  temporal  and  spatial  scales.  Even 


if  we  could,  uncertainty  in  boundary 
conditions  and  forcing  would  require  us 
to  consider  even  the  most  deterministic 
model  results  in  a  statistical  sense.  One 
response  to  this  dilemma  is  to  adopt  a  hi¬ 
erarchy  of  models  that  inform  each  other 
across  a  range  of  domains  and  scales.  The 
Adriatic  modeling  work  presented  here 
is  an  example  of  this  approach.  We  use 
multiple  nested  meteorological  models  to 
transfer  information  from  global  scales 
to  local  scales  (consisting  of  a  few  square 
kilometers).  The  coupling  of  atmospher¬ 
ic,  hydrologic,  wave,  and  oceanographic 
models  allows  us  to  evaluate  the  impor¬ 
tance  of  regional  and  local  forcing  on  the 
detailed  physics  of  the  bottom  boundary 


layer,  which  we  can  examine  in  detail 
with  the  gravity-flow  models.  Coupled 
circulation  and  sediment-transport  mod¬ 
els  allow  us  to  verify  processes  against 
measurements  during  specific  events 
and  over  the  course  of  a  winter  transport 
season.  Stratigraphic  models  let  us  ex¬ 
tend  the  effect  of  annual  processes  over 
decades  and  centuries  and  factor  in  long¬ 
term  changes  in  sea  level  and  sediment 
supply.  In  the  end,  we  hope  to  achieve  a 
cascade  of  models  that,  although  they  are 
calibrated  at  small  scales  and  over  short 
periods,  still  allow  us  to  build  insight 
about  cumulative  effects  because  they  are 
based  on  a  fundamentally  correct  repre¬ 
sentation  of  the  processes. 
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OBSERVATIONS 

EuroSTRATAFORM  researchers  de¬ 
ployed  instrumented  tripods  and  moor¬ 
ings  for  nine  months  (November  2002 
to  May  2003);  made  water-column 
measurements  in  the  western  Adriatic 
coastal  current  (WACC)  during  voyages 
in  November,  February,  and  May;  made 
short-term  (~2  day)  studies  of  particle 
dynamics;  and  analyzed  cores  to  de¬ 
termine  characteristics  of  near-surface 
sediments  and  accumulation  rates  over 
time  scales  ranging  from  weeks  to  mil¬ 
lennia.  Participants  in  the  ACE,  EACE, 
and  WISE  programs  deployed  acoustic 
Doppler  current  meter  moorings  in  the 
northern  Adriatic  and  collected  hydro- 
graphic  data  in  September  to  October 
2002  and  late  May  2003  (Figure  1). 
DOLCE  VITA  scientists  tracked  scores 
of  drifters  and  measured  water-column 
sections  of  temperature,  salinity,  fluo¬ 
rescence,  and  light  attenuation  during 
two  cruises  in  February  and  May/June. 
ADRICOSM  researchers  conducted 
monthly  surveys  of  water-column  prop¬ 
erties  in  several  key  near-coastal  regions. 
These  measurements  provide  a  rich  set 
of  data  for  initializing,  forcing,  and  criti¬ 
cally  evaluating  numerical  models  of  the 
major  processes  that  ultimately  write  the 
stratigraphic  record  in  the  Adriatic  Sea. 

MODELING 

We  have  linked  a  sequence  of  models  that 
(1)  characterize  river  flow  and  sediment 
discharge,  (2)  estimate  time-varying  me¬ 
teorological  conditions  and  waves,  (3) 
compute  circulation  and  sediment  trans¬ 
port  for  a  nine-month  winter  period,  (4) 
simulate  the  remobilization,  transport, 
and  deposition  of  sediments  during  fluid 
mud  events,  and  (5)  simulate  the  devel¬ 


opment  of  stratigraphy. 

Winter  circulation  in  the  Adriatic  Sea 
is  strongly  controlled  by  local  wind  and 
river  input  (e.g.,  Orlic  et  al.,  1994;  Cush- 
man-Roisin  et  al„  2001),  so  it  is  an  ex¬ 
cellent  natural  laboratory  for  numerical 
modeling  studies.  Time  series  of  the  sup¬ 
ply  of  freshwater  and  sediment  from  the 
Po  River  and  smaller  Apennine  rivers  are 
required,  as  are  fields  of  meteorological 
conditions  (such  as  winds,  temperature, 
humidity,  precipitation)  that  control 
air-sea  exchanges  of  momentum,  heat, 
and  freshwater.  Waves  play  a  critical  role 
in  sediment  resuspension,  so  accurate 
knowledge  of  near-bottom  wave  motions 
is  necessary.  Once  sediment  is  resuspend¬ 
ed  by  waves  or  delivered  by  rivers,  region¬ 
al  currents  transport  it;  hence,  a  detailed 
description  of  circulation  is  required, 
particularly  during  transport  events. 

The  data  collected  in  2002  to  2003  and 
the  suite  of  models  now  available  have 
allowed  us  to  estimate  depositional  pat¬ 
terns  and  compare  them  with  the  long¬ 
term  accumulation  patterns  observed  in 
sediment  cores  and  seismic  data. 

Meteorology 

Meteorological  forcing  for  the  wave  and 
circulation  models  was  output  from  two 
models:  LAMI  (Limited  Area  Model 
Italy)  and  COAMPS™  (Coupled  Ocean- 
Atmosphere  Mesoscale  Prediction  Sys¬ 
tem).  Both  solve  finite-difference  ap¬ 
proximations  of  the  three-dimensional, 
non-hydrostatic  equations  that  govern 
atmospheric  motions.  LAMI  is  an  Ital¬ 
ian  operational  model  for  medium-  and 
small-scale  weather  prediction  based  on 
a  model  developed  by  the  German  me¬ 
teorological  service  (Deutscher  Wetter- 
dienst)  (Steppler  et  al.,  2003).  COAMPS 


(Hodur  and  Doyle,  1999)  is  the  opera¬ 
tional  mesoscale  prediction  system  used 
in  support  of  the  U.S.  Navy,  and  it  can 
be  coupled  with  the  Navy  Coastal  Ocean 
Model  (NCOM).  Both  atmospheric 
models  were  nested  into  regional  and 
global  weather  models.  COAMPS  model 
results  were  generated  more  frequently 
and  on  a  finer  grid  ( 1  hr  at  4  km,  com¬ 
pared  with  3  hrs  at  7  km  for  LAMI),  and 
COAMPS  runs  assimilated  data  from 
radiosondes,  surface  satellite,  and  aircraft 
observations.  In  “one-way”  coupled  runs, 
the  atmospheric  models  provided  forc¬ 
ing  for  the  ocean  models  and  measured 
or  modeled  sea-surface  temperature  was 
used  to  estimate  air-sea  heat  exchange. 
When  “two-way”  coupled  with  NCOM, 
COAMPS  forced  the  ocean  model  with 
momentum  and  heat  flux  and  NCOM 
supplied  sea-surface  temperatures  as 
feedback  to  COAMPS. 

River  Discharge 

Rivers  are  a  primary  source  of  sediment, 
and  the  freshwater  they  deliver  influences 
coastal  circulation;  thus,  correct  repre¬ 
sentation  of  river  contributions  in  the 
model  is  vital.  We  set  flow  according  to 
daily-average  measurements  for  the  Po, 
Pescara,  and  Biferno  Rivers.  We  set  the 
flow  of  other  rivers  to  monthly  mean  val¬ 
ues  using  climatological  estimates.  Sedi¬ 
ment  discharge  was  estimated  using  the 
one-dimensional  flow  model  HydroTrend 
(Syvitski  et  al.,  1998),  which  creates  syn¬ 
thetic  river  discharge  and  sediment-load 
time  series  as  a  function  of  climate  trends 
and  basin  morphology.  The  sediment 
rating  curves  generated  by  HydroTrend 
were  combined  with  measured  discharges 
to  generate  time  series  of  sediment  con¬ 
centrations  in  the  Po  River  for  input  to 
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2002-2003  model  simulations.  River-sup- 
plied  sediment  was  broadly  characterized 
as  either  silt,  with  a  low  settling  velocity 
of  0.1  mm/s,  or  flocculated  finer  material, 
with  an  aggregate  settling  velocity  of  1 
mm/s  (Fox  et  al.,  2004;  in  press). 

Waves 

Waves  were  simulated  using  SWAN 
(Simulating  WAves  in  the  Nearshore), 
which  models  the  evolution  of  wave  en¬ 
ergy  imparted  by  wind  as  waves  propa¬ 
gate  across  the  sea;  are  transformed  (in 
frequency  and  direction)  by  shoaling, 
refraction,  and  wave-wave  interactions; 
and  are  dissipated  by  whitecapping  and 
bottom  friction  (Booij  et  al.,  1999).  We 
modeled  waves  in  the  Adriatic  Sea  on  a 
rectangular  grid  with  ~2-km  spacing. 
Evolution  of  the  non-stationary  wave 
field  was  calculated  at  15-minute  time 
steps,  driven  by  wind  fields  provided  at 
three-hour  intervals  and  ~7-km  spacing 
from  the  LAMI  model  output. 

Circulation  and  Suspended- 
Sediment  Transport 
Oceanic  circulation  and  sediment  trans¬ 
port  were  simulated  using  ROMS  (Re¬ 
gional  Ocean  Modeling  System)  (more 
information  available  at  http://marine. 
rutgers.edu/po/models/roms/) .  ROMS 
solves  finite-difference  approximations 
of  the  three-dimensional  Reynolds-aver¬ 
aged  equations  for  conservation  of  mass, 
momentum,  and  heat  using  a  two-equa¬ 
tion  submodel  for  turbulent  mixing. 
Effects  of  waves  in  the  bottom  bound¬ 
ary  layer,  which  include  increased  drag 
and  bottom  stresses,  were  parameterized 
with  a  modified  version  of  the  Grant  and 
Madsen  (1979)  model.  Sediment  sup¬ 
plied  by  rivers  or  resuspended  from  the 


seafloor  was  advected,  vertically  mixed, 
and  allowed  to  settle  in  the  model,  which 
also  records  changes  in  upper  stratigra¬ 
phy  (to  ~0.1  m)  of  the  seafloor  caused  by 
erosion  or  deposition.  Wind-driven  cir¬ 
culation,  mixing,  and  heating  or  cooling 
of  surface  waters  were  calculated  using 
bulk  flux  algorithms  with  input  from  the 
atmospheric  models.  Initial  water  tem¬ 
peratures  and  salinities  were  set  by  in¬ 
terpolating  from  hydrographic  data  col¬ 
lected  during  September  2002  surveys. 

A  range  of  model  runs  were  performed 
on  a  ~4-km  grid  with  20  vertical  levels 
using  both  atmospheric  models  and  al¬ 
ternative  parameterizations  for  various 
processes  to  explore  the  fate  of  sediment 
supplied  by  rivers  or  resuspended  from 
existing  deposits. 

Gravity  Flows  of  Fluid  Mud 
Gravity-flow  processes  on  the  Po  River 
delta  front  were  simulated  with  a  one-di¬ 
mensional  (cross-shelf)  time- dependent 
model  that  resolves  the  suspended  sedi¬ 
ment  concentration  in  the  thin  wave- 
current  boundary  layer  and  calculates 
transport  and  deposition  associated 
with  gravity-driven  downslope  flows. 

The  model  was  forced  with  observed 
waves  and  currents,  and  assessed  for  skill 
against  near-bottom  measurements  of 
flow  and  suspended-sediment.  A  two-di¬ 
mensional  model  (across  shelf  and  along 
shelf)  based  on  several  simplifying  as¬ 
sumptions  was  also  used  to  predict  the 
spatial  distribution  of  the  flood  deposit 
on  the  Po  delta  (see  Nittrouer  et  al.,  this 
issue).  These  models  allow  us  to  evalu¬ 
ate  the  displacement  of  flood  deposits 
from  their  initial  depocenters  beneath 
the  river  plume  to  their  eventual  locus  of 
long-term  accumulation. 


Stratigraphy 

Stratigraphy  was  considered  on  two 
scales:  detailed  (upper  2  to  10  cm)  and 
shallow  (upper  ~2  to  30  m).  Two  mod¬ 
els,  HSM-3D  and  CST  (Coastal  Sedi¬ 
ment  Tract),  have  been  developed  to 
allow  diagnosis  and  predictions  at  these 
scales.  HSM-3D  is  a  three-dimensional, 
time-dependent  numerical  model  for 
hydrodynamics,  sediment  dynamics,  and 
morphodynamics.  It  includes  suspended 
transport,  bedload,  and  the  effect  of  grav¬ 
ity  flows  in  the  bottom  boundary  layer 
and  incorporates  evolving  bottom  mor¬ 
phology  caused  by  erosion  or  deposition. 
The  model  has  variable  grid  spacing  in 
the  water  column  and  fixed  bins  for  stra¬ 
tigraphy,  which  can  extend  several  meters 
below  the  seafloor.  We  embedded  the 
HSM-3D  grid  in  the  ROMS  domain  to 
resolve  length  scales  of  ~0.1  km,  while  in¬ 
corporating  information  about  the  larger 
scale  mean  circulation  and  wave  fields. 

The  CST  model  simulates  sedimen¬ 
tation  over  decades  and  millennia  by 
parameterizing  sediment  dynamics  in 
terms  of  time-averaged  and  spatially  av¬ 
eraged  relationships  that  are  consistent 
with  long-term  system  forcing  (Cowell  et 
al.,  in  press).  The  CST  model  represents 
a  number  of  sedimentary  systems  in¬ 
cluding  the  river  watershed,  estuary,  in¬ 
let,  beach/surf  zone,  and  shoreface/shelf 
systems  (see  Pratson  et  al.,  this  issue). 

RESULTS 

Meteorology  and  Wave  Models 
Wind  fields  and  air-sea  heat  fluxes  de¬ 
rived  from  both  meteorological  models 
COAMPS  and  LAMI  matched  each  other 
and  field  measurements  (e.g„  winds 
measurements  offshore  Venice)  in  a 
broad  sense  (correlation  ~0.65;  error  in 


Oceanography  I  Vol.17,  No.4,  Dec.  2004  63 


amplitude  response  ~10  percent), but 
differed  in  the  details  (Cavaleri,  1999). 
Compared  to  the  coarser  resolution  runs, 
the  highest-resolution  COAMPS  nest  (4 
km)  demonstrated  better  skill  in  produc¬ 
ing  the  spatial  pattern  of  Bora  winds, 
and  this  led  to  improvements  in  mod¬ 
eled  wind-forced  ocean  currents  (Pullen 
et  al.,  2003).  These  wind  patterns  also 
seem  to  be  captured  by  LAMI.  Some  of 
the  small-scale  structure  is  deterministic 
and  reflects  improved  resolution  of  oro¬ 
graphic  effects,  but  some  is  random  and 
can  degrade  single-point  correlations  be¬ 
tween  model  results  and  measurements 
(Signell  et  al.,  in  press).  Wave  fields  and 
ocean  heat  content  naturally  integrate 
meteorological  forcing;  good  agreement 
among  measured  and  modeled  ocean 
temperatures,  waves,  and  currents  in  the 
Adriatic  Sea  model  suggests  the  meteo¬ 
rological  models  were  generally  accurate. 

COAMPS/NCOM  simulations  dem¬ 
onstrate  the  advantages  of  two-way  cou¬ 
pling,  where  more  accurate  water  tem¬ 
peratures  improve  the  parameterizations 
for  atmospheric  boundary  layer  struc¬ 
ture  and  heat  fluxes  (Figure  3).  But  the 
one-way  LAMI/ROMS  simulations  seem 
to  capture  key  aspects  of  the  oceanic 
response  to  atmospheric  forcing,  and 
produced  heat  fluxes  and  sea-surface 
temperatures  that  fell  within  the  uncer¬ 
tainty  of  preliminary  estimates  derived 
from  five  cloud-free  satellite  images  and 
in  situ  temperature  data  (Figure  3b).  At¬ 
mospheric  forcing  is  the  most  important 
input  to  our  sediment-transport  models 
in  the  Adriatic  Sea  and  an  ongoing  re¬ 
search  topic. 

The  wave  climate  in  winter  2002  to 
2003  was  dominated  by  a  few  Sciroccos 
and  a  series  of  Bora  wind  events.  Maxi- 
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Figure  3.  (a)  As  cold  Bora  winds  cross  the  Adriatic,  they  fiow  over  surface  wa¬ 
ters  that  are  ~10  degrees  warmer.  This  induces  vigorous  air-sea  exchange  as 
the  atmosphere  extracts  heat  from  the  ocean.  The  vertical  contour  plot  shows 
COAMPS/NCOM  two-way  coupled  potential  temperature  at  the  beginning  of 
a  Bora,  (b)  Bulk  heat  transfer  rates  calculated  by  ROMS  fall  within  error  bounds 
around  fluxes  estimated  from  satellite  and  in  situ  temperature  data,  but  LAMI 
algorithms  overpredict  rate  of  heat  transfer. 
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mum  modeled  waves  reached  5  m  (sig¬ 
nificant  wave  height)  in  open  waters  of 
the  northeast  Adriatic  during  a  Scirocco 
that  inundated  Venice  lagoon  with  the 
sixth  highest  aqua  alta  (storm  surge) 


in  the  last  60  years  (more  information 
available  at  Centro  Previsioni  e  segnala- 
zioni  Maree,  Comune  di  Venezia;  http:// 
www.comune.venezia.it/maree/).  Bora 
events  were  remarkable  for  their  spa¬ 


tial  structure  (Figure  4).  Comparisons 
of  SWAN  model  hindcasts  with  Italian 
buoy  records  indicate  the  model  cap¬ 
tures  most  of  the  features  in  the  records 
of  wave  height.  Modeled  near-bottom 


wave  orbital  velocities  also  agree  fairly 
well  with  PASTA  tripod  measurements, 
but  there  are  some  intriguing  differ¬ 
ences  during  some  events  (Figure  4)  that 
may  be  related  to  incorrect  wind  forcing 
or  problems  with  swell  dissipation  in 
SWAN  (see  Rodgers  et  al.,  2003). 


Circulation,  Sediment-Transport, 
and  Stratigraphy  Models 
Modeled  circulation  in  the  Adriatic 
agreed  well  with  the  extensive  data  col¬ 
lected  in  2002  to  2003.  The  mean  flow, 
driven  by  winds  and  buoyancy  input, 
matches  circulation  described  by  ear¬ 


lier  drifter  studies  (Poulain,  2001)  and 
recent  modeling  efforts  (e.g.,  Pullen  et 
al„  2003).  Detailed  comparison  during 
Bora  events  (Figure  5)  indicates  that  the 
model  captures  complex  circulation  pat¬ 
terns.  Modeled  transport  in  the  WACC 
was  compared  with  estimates  made  from 


Wind  at  10  m  (m/s)  Currents  at  -  5  m  (m/s) 


Figure  5.  (a)  Bora  winds  form  a  distinct 
pattern  across  the  northern  Adriatic 
as  they  funnel  through  topography  in 
Croatia,  (b)  The  resulting  circulation  is 
complex,  as  shown  here  by  simulated 
mean  flow  over  four  days  at  water  depth 
of  5  m.  Model  results  from  ROMS  gener¬ 
ally  correspond  well  with  measurements 
from  the  ACE  array  of  acoustic  Dop¬ 
pler  current  meters  (magenta  arrows), 
(c)  Drifter  tracks  simulated  with  ROMS 
using  wave-enhanced  near-surface  mix¬ 
ing  (solid  green)  match  tracks  measured 
during  DOLCE  VITA  (black)  better  than 
tracks  simulated  without  wave-en¬ 
hanced  mixing  (dashed  green). 
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Figure  6.  The  western  Adriatic  current  causes  southward  water  transport  along  the  Italian  coast. 
Model  results  (blue)  match  the  timing  and  often  the  magnitude  of  transport  estimated  from 
an  array  of  four  acoustic  Doppler  current  profiles  on  a  transect  near  Senigallia  (see  Figure  5b). 


four  acoustic  Doppler  current  profiles 
(ADCPs)  off  Senigallia  (Figure  6).  The 
timing  of  modeled  transport  events 
agreed  closely  with  observations,  but  the 
peaks  of  the  events  were  underestimated. 
Overall,  the  modeled  mean  transport 
was  about  20%  lower  than  the  mean  in¬ 
ferred  from  measurements  on  this  tran¬ 
sect.  Sediment  transport  is  strongly  non¬ 
linear,  however,  and  the  impact  of  errors 
in  circulation  during  specific  events  is 
under  investigation. 

Model  simulations  of  river-supplied 
sediment  naturally  indicated  that  rapid¬ 
settling  material  tends  to  accumulate 
near  the  river  mouths  (particularly  the 
Po),  and  that  unflocculated  material  is 
transported  farther  along  the  coast  and 
accumulates  in  elongated  deposits  near 
the  shelf  edge  (Figure  7).  The  depo- 
sitional  pattern  produced  by  a  model 
simulation  for  winter  2002  to  2003  mim¬ 
ics  aspects  of  the  long-term  depositional 
pattern  observed  in  Holocene  stratigra¬ 
phy.  Specifically,  deposition  on  the  Po 
delta,  near  the  Apennine  river  mouths, 
and  near  the  shelf  edge  north  of  the 
Gargano  promontory  occurs  where  late 
Holocene  deposits  are  thickest.  Typical 
results  from  the  HSM-3D  model  showed 
depositional  patterns  that  indicate 
alongshore  transport  produces  coalesc¬ 
ing  clinoforms  off  the  Apennine  rivers. 
Larger-scale  features  that  develop  over 
millennia  are  predicted  with  the  CST 
model  (Pratson  et  ah,  this  issue).  The 
encouraging  agreement  suggests  that  the 
linked  system  of  models  incorporate  the 
key  mechanisms  for  transport  and  deliv¬ 
ery  of  sediment  from  shelf  environments 
to  depocenters  located  farther  south  in 
deeper  water.  Work  is  ongoing  to  inves¬ 
tigate  the  sensitivity  of  these  patterns 


to  winds,  settling  velocity,  variations  in 
erosion  formulae,  and  other  parameters. 
These  sensitivity  studies  will  help  us  un¬ 
derstand  the  weak  links  in  the  overall 
modeling  system. 

SEDIMENT-TRANSPORT 

MECHANISMS 

Measurements  and  model  simulations  of 
winter  conditions  in  the  Adriatic  showed 
that  near  the  mouth  of  the  Po  River,  rap¬ 
id  accumulation  after  floods  produces 
ephemeral  deposits  of  sediment  that  are 
subsequently  remobilized  by  waves  to 
form  density  flows.  This  was  observed 
on  three  occasions  in  2002  to  2003,  and 
these  fluid-mud  flows  moved  sediment 
from  ephemeral  shallow-water  deposits 
beneath  the  river  plume  to  deeper  water. 
These  results  lead  us  to  believe  that  fluid 
mud  processes  may  control  the  geometry 
of  deltaic  deposits;  EuroSTRATAFORM 
researchers  are  developing  models  to  cal¬ 
culate  the  resulting  depth  profiles. 

The  mechanisms  that  move  Po  River 


sediment  hundreds  of  kilometers  to  de¬ 
pocenters  off  the  Gargano  promontory 
must  be  associated  with  a  larger-scale 
phenomena  like  the  western  Adriatic 
coastal  current.  The  WACC  is  partially 
buoyancy  driven,  a  forcing  that  increases 
during  floods  when  suspended-sedi¬ 
ment  supply  is  highest.  Observations 
and  models  demonstrate  that  Po  water 
flows  around  the  Gargano  Promontory, 
but  because  the  sediment  aggregates  and 
settles  rapidly,  it  is  unlikely  that  much  re¬ 
mains  in  suspension  for  the  entire  jour¬ 
ney.  It  appears  that  Po  sediments  must 
travel  southward  in  a  series  of  episodic 
transport  events,  and  that  wave-induced 
resuspension  is  important  in  these  hops 
(Wang  and  Pinardi,  2002). 

Bora  winds  generate  large  waves  in  the 
western  Adriatic  and  enhance  flow  in  the 
WACC,  so  it  is  reasonable  to  expect  that 
correlated  wave  resuspension  and  stron- 
ger-than-average  southward  flow  in  the 
WACC  combine  to  generate  southward 
sediment  flux.  However,  our  data  and 
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Figure  7.  (a)  Results  of  sediment-transport  simulations  for  nine  months  (November  2002  -  May  2003)  reveal  patterns  of  deposition 
for  (a)  fast-settling  aggregates  (1  mm/s)  and  (b)  slow-settling  particles  (0.1  mm/s)  supplied  by  rivers.  The  depositional  patterns  are 
dominated  by  local  deposition  of  material  from  the  Po  River  in  the  northern  Adriatic  and  by  southward  dispersal  of  Apennine  sedi¬ 
ment  along  much  of  the  Italian  coast.  Deposition  is  enhanced  north  of  Ancona  and  the  Cargano  promontory,  similar  to  the  pattern 
seen  in  the  Holocene  thicknesses. 


model  results  show  only  a  weak  correla¬ 
tion  between  wave  height  and  southward 
flow  velocity  in  the  WACC  on  the  Chien- 
ti  transect.  The  best  correlation  between 
WACC  flow  at  Chienti  and  modeled 
wind  is  with  wind  near  Trieste,  nearly 
200  km  from  the  Chienti  transect.  This 
is  evidence  that  wind-driven  circulation 
in  the  northern  Adriatic  has  a  very  large 
scale,  and  that  waves  are  not  necessarily 
collocated  with  Bora-enhanced  south¬ 
ward  flow.  The  model  indicates  that 
Scirocco  winds  can  generate  large  waves 


that  resuspend  sediment  without  revers¬ 
ing  flow  in  the  WACC,  suggesting  that 
southward  transport  is  likely  to  prevail 
regardless  of  the  timing  or  mechanisms 
for  resuspension. 

The  system  of  models  we  have  ap¬ 
plied  to  the  Adriatic  Sea  incorporates 
processes  at  scales  ranging  from  hours 
and  kilometers  up  to  millennia  and  ba¬ 
sin  dimensions.  Large  collaborative  field 
efforts  such  as  those  in  the  Adriatic  Sea 
during  2002  to  2003  are  the  only  viable 
way  to  obtain  enough  synoptic  data  on 


waves,  currents,  and  sediment  to  evalu¬ 
ate  such  complex  models.  Although 
research  is  ongoing  to  assess  the  skill  of 
these  models,  results  from  the  simula¬ 
tions  are  already  helping  us  gain  a  better 
understanding  of  the  mechanisms  that 
form  deposits  on  the  margins  of  epicon¬ 
tinental  seas. 
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